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We describe experiments and simulations demonstrating the propulsion of a neutrally-buoyant 
swimmer that consists of a pair of spheres attached by a spring, immersed in a vibrating fluid. The 
vibration of the fluid induces relative motion of the spheres which, for sufficiently large amplitudes, 
can lead to motion of the center of mass of the two spheres. We find that the swimming speed 
obtained from both experiment and simulation agree and collapse onto a single curve if plotted as a 
function of the streaming Reynolds number, suggesting that the propulsion is related to streaming 
flows. There appears to be a critical onset value of the streaming Reynolds number for swimming 
to occur. We observe a change in the streaming flows as the Reynolds number increases, from that 
generated by two independent oscillating spheres to a collective flow pattern around the swimmer 
as a whole. The mechanism for swimming is traced to a strengthening of a jet of fluid in the wake 
of the swimmer. 


The mechanism by which self-propulsion through a 
fluid is achieved has fascinated scientists of many dis¬ 
ciplines, and the public alike, for aesthetic, practical 
and fundamental scientific reasons nHa. In biology and 
biomechanics the mechanisms behind the way organisms 
swim gives insight into their biological function and pur¬ 
pose [Ilia EE]. Recently, the design of efficient “robots” 
able to navigate themselves through various fluids has be¬ 
come an important technological and medical challenge 
that brings together elements of physics, chemistry, biol¬ 
ogy, engineering and fluid mechanics iHlo]. Microscopic 
artificial swimmers have been proposed for use in tar¬ 
geted drug-delivery, see for example [IIHI3]. 

Purcell’s scallop theorem states that at zero Reynolds 
number an object cannot swim using a time-reversible 
stroke: it will end up going back and forth with no net 
displacement [3]. Many types of small creatures, for ex¬ 
ample, insects and aquatic invertebrates, swim at inter¬ 
mediate Reynolds numbers (1-100) [M]. In these cases, 
time-reversal symmetry is broken by non-linearities in the 
fluid dynamics rather than by the nature of the stroke. 
For such swimmers, an interesting question arises: how 
does the motion evolve as the Reynolds number is in¬ 
creased from zero? It has been argued that symmet¬ 
rical objects with symmetrical strokes such as flapping 
wings have an onset for motion at a critical Reynolds 
number HMH, whereas asymmetrical objects or strokes 
have a continuous transition as the Reynolds number is 
increased m- 

A central problem when designing a practical arti¬ 
ficial swimmer is how to get energy into the system. 
Methods based on electromagnetic or chemical actuation 
have been developed m and currently there is inter¬ 
est in using acoustic techniques to generate propulsion 
through the oscillation of entrapped air bubbles [20l [2T] . 
Vladimirov proposed an alternative mechanism that may 


lead to swimming based on a deformable object which is 
neutrally buoyant, but composed of coupled spheres with 
different sizes and densities [22]. Such an object can gen¬ 
erate relative motion of its parts if immersed in a vibrat¬ 
ing fluid; this motion may lead to swimming. However, 
his calculations in the absence of fluid and particle inertia 
predicted that such an object will not swim if subjected 
to unidirectional oscillation. Here we pose the question: 
can an experimental realisation of this object be made to 
swim at higher Reynolds numbers, and, if so, what is the 
method of propulsion and the nature of the transition to 
swimming? 

In this Letter we describe experiments and simulations 
demonstrating the propulsion of a pair of spheres at¬ 
tached by a spring, immersed in a vertically vibrating 
fluid. We consider two particular realisations of this ob¬ 
ject: one with unequal-sized spheres and the other with 
equal-sized spheres. In both cases, the density of the 
spheres is different from one another and from the liq¬ 
uid in which they are immersed, however, the object 
as a whole is neutrally-buoyant. We find that both de¬ 
signs swim for sufficiently high amplitudes of vibration; 
the unequal-sized spheres swim upwards, in the direc¬ 
tion of the larger, less dense sphere, whereas the equal¬ 
sized spheres swim downwards, in the direction of the 
higher density sphere. The data for the swimming speed 
are found to collapse both in experiment and simulation 
when scaled appropriately with the streaming Reynolds 
number, suggesting that the streaming flows induced by 
fluid non-linearities play a central role [23|. Furthermore, 
the apparent onset of motion appears to be governed by 
a critical value of the streaming Reynolds number. The 
mechanism for propulsion is traced to a change in the 
topology of the streaming flows that transition from those 
of two noninteracting spheres when the dimer is station¬ 
ary, to a collective flow around the object, at the apparent 
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FIG. 1. Main panel shows the experimental data collapse 
for the unequal-sized swimmer which swims upwards. The 
driving conditions are the following: blue stars / = 65Hz, 
light blue triangle down / = 75Hz, pink hexagons / = 85Hz, 
cyan circles / = 95Hz, yellow triangle up / = 105Hz, green 
triangles right / = 125Hz, green diamond / = 135Hz. In 
all cases the viscosity was 1.2mm^/s except for one data 
set (turquoise diamond / = 135Hz) where the viscosity was 
2.5mm^/s. Simulations for F between 2 — 20 and frequencies 
/ = 65, 75,100,125Hz are shown by the red plus symbols for 
comparison. The lower inset shows a close-up of the experi¬ 
mental data (/ = 125Hz) shown in the main panel near the 
apparent onset. The upper inset shows a photograph of the 
swimmer when stationary. 


onset of swimming. The flow field shows a strengthening 
of a jet of fluid behind the swimmer. 

The dimers were constructed from two spheres joined 
together by a small coil of wire. Examples of the asym¬ 
metric and symmetric dimers are shown in the insets to 
Figs. 1 and 2 respectively. Details of their construction 
and the experimental set-up are given in Supplemental 
Material m- The dimers were designed so that they 
could be made neutrally buoyant in a salt-water solution. 
The solution was vibrated vertically at a given frequency, 
/, amplitude, A. The dimensionless acceleration of the 
cell F = A{27rf)‘^/g was varied between 1 — 20, where g 
is the gravitational acceleration. The frequency ranged 
from 30Hz to 135Hz. 

As the cell vibrated, each sphere had a different am¬ 
plitude and phase relative to the fluid motion due to dif¬ 
ferences in the size and/or densities of the two spheres. 
At low amplitudes of vibration of the cell the spheres 
oscillated vertically, but no net time-averaged motion 
of the center of mass of the spheres could be observed 
within experimental error. Beyond a certain thresh¬ 
old the dimer started to move; increasing the amplitude 
made the dimer swim faster. 

To obtain the velocity of the dimer, the vibration was 
initiated abruptly under fixed F and /, and the motion 
of the dimer was filmed using a high-speed camera. The 


movies show that the separation of the spheres varied si¬ 
nusoidally (indicating that the coiled wire acted as a lin¬ 
ear spring to a good approximation). From such movies 
the steady-state velocity of the dimer, and the relative 
amplitude of the two spheres with respect to each other, 
Ar^ was obtained. Ar is the amplitude of the relative 
motion of the two spheres that comprise the swimmer. 
Note that A^ and the driving amplitude of the cell are 
different; A^ increases approximately linearly with A. As 
far as the motion of the spheres is concerned, in the rest 
frame of the cell, Ar and / are the only relevant driving 
parameters. As can be seen from the movies [24], the mo¬ 
tion of the spheres was predominantly along the vertical 
line through their centers; there was very little sideways 
‘waggling’ movement. 

Fig. 1 shows the data obtained for the two asymmetric 
dimers, which swim upwards in the direction of the larger 
sphere. The data collapse (within the scatter) when plot¬ 
ted in terms of the dimensionless combinations v/fL^ and 
the streaming Reynolds number Rcs = Here L 

is the diameter of the larger sphere and 6 = {y/2'Kf)^ 
is the viscous length in terms of the kinematic viscosity 
V. Rcs is one of three dimensionless ratios that can be 
defined from the length scales L and d and charac¬ 
terises the time-averaged (steady) flow [23 HI] • In onr 
experiments L ^ S which results in a configuration of 
the time-averaged flow around each sphere that has in¬ 
ner and outer loops PH EO]. The data are consistent 
between the measurements obtained from two nominally 
identical, asymmetric dimers, indicating that small dif¬ 
ferences in construction such as variations in the shape 
of the loop of wire and of the shape and amount of glue 
have little effect. The collapse in terms of Rcs shows 
that the motion is related to streaming flows generated 
by the vibration of the dimer. The lower inset shows 
data taken at low amplitudes of vibration and suggests a 
sharp increase in velocity at Res ~ 20. 

Fig. 2 shows the behaviour of the equal-sized dimer. In 
this case it moves in the opposite direction, downwards, 
with the heavier sphere at the front. The data illustrates 
that the speed and direction of motion depends on the 
densities and sizes of the spheres, as well as the gap be¬ 
tween them; if the two spheres are sufficiently far apart, 
the dimer will not swim. The main reason for consider¬ 
ing this system is that it is arguably one of the simplest 
objects that can be made to swim. Note that it was 
difficult to design dimers made of equal-sized spheres of 
different densities that can be made neutrally buoyant in 
salt solutions, and have sufficient mass difference between 
the spheres to generate enough relative motion to induce 
swimming. Hence the relative paucity of data compared 
to that obtained for the unequal-sized dimers. 

We now ask the question: what is the cause of the mo¬ 
tion? To address this we first imaged the flow using tracer 
particles illuminated by a planar laser sheet in the plane 
containing the centres of the two spheres of the dimer. 
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FIG. 2. Main panel shows the experimental data collapse 
for the equal-sized swimmer. The driving conditions are the 
following: blue stars / = 30Hz and green triangles right / = 
35Hz. Simulations for F between 2—12 and frequency / = 
30Hz are shown by the red plus symbols for comparison. The 
inset shows a photograph of the swimmer when stationary. 


Fig. 3 (a) shows a photograph of the asymmetric dimer 
taken with an exposure time of one period of oscillation, 
revealing the motion of the tracer particles. In this image 
the dimer is close to the onset of motion. A downward jet 
originates from the vicinity of the lower sphere [24] . Sim¬ 
ilar behaviour was found for equal-sized spheres, except 
that the strong jet was generated by the upper, lighter 
sphere, causing the swimmer to swim downwards. 

In order to investigate the motion of the spheres and 
the fluid in detail we used simulations which were based 
on an embedded boundary method described previ¬ 
ously [3QH34] . The fluid was assumed to obey the Navier- 
Stokes equations which were discretised on a staggered 
mesh [35] and solved using the projection method [36] 
to ensure incompressibility of the fluid. The interac¬ 
tion between the fluid and the rigid spheres was achieved 
through the template model, which introduces a two-way 
coupling between the particles and the fluid [33]. The 
spheres were joined by a linear spring as in the experi¬ 
ments. An equal and opposite force was applied vertically 
to the spheres to mimic the effects of static buoyancy, 
rather than imposing the effect of gravity directly on the 
fluid. The influence of vibration was introduced by ap¬ 
plying a sinusoidal acceleration to the fluid and particles, 
so that the simulations were carried out in the frame of 
reference of the vibrated cell. 

The computational parameters of the swimmer (size, 
density and gap) and fluid (viscosity and density) were 
chosen to match the experiments. However, any interac¬ 
tion of the wire with the fluid was ignored and the dimers 
were assumed to be made of perfect spheres. Details of 
the parameters used are given in the Supplemental Mate¬ 
rial. One difference between experiment and simulation 



FIG. 3. Illustrations of the fluid flows generated by the vi¬ 
bration of the spheres from experiment and simulation. Panel 
(a) shows an image taken from experiment showing the flow 
around the spheres. The arrow illustrates the direction of a 
jet of fluid evident from the movies (supplemental informa¬ 
tion [13]). Panels (b) and (c) show the direction of the time- 
averaged velocity field (i.e. the normalised velocity vectors) 
in the plane of the spheres. In (b) the swimmer is stationary 
(Res = 15) while in panel (c) it is swimming (Rcg = 60). 
These figures illustrate the change in topology of the flows as 
the amplitude of vibration increases. Note that the magni¬ 
tude of the flows is much greater around the smaller sphere 
than around the larger sphere, as seen in panel (a). 


is that the simulated cells are smaller due to computa¬ 
tional limitations. Examples of the simulated data are 
shown in Figs. 1 and 2 by the large red plus symbols. 
There is clearly good agreement between the simulations 
and experiment despite the numerical limitations arising 
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FIG. 4. Data collapse from simulations confirming the scaling 
behaviour for different viscosities (red 1.2 x 10 ®m^ /s, blue 
2 X 10“®m^/s, green 3 x 10“®m^/s). Each data set includes 
simulations for F between 2 — 20 and / = 65, 75,100,125Hz. 
The lower inset shows typical trajectories of a swimmer that 
stays stationary for F = 2, f — 65Hz (red line) and one that 
swims for F = 12, f — 75Hz (blue line). The upper left inset 
is a snapshot from simulations showing the swimmer and the 
simulated cell. 

from the simulated cell size and possible fluid lattice ef¬ 
fects. 

The simulations allow us to determine in more detail 
the fluid flows generated by the motion of the spheres in¬ 
duced by the vibration of the cell. This flow is best illus¬ 
trated by plotting the direction of the velocity field in the 
vertical plane through the center of the two spheres. Ex¬ 
amples of these flows for the two unequal-sized spheres, 
time-averaged over a cycle, are shown in Fig. 3 (b) and 
(c). At low amplitudes. Fig. 3 (b), there are two outer 
vortex rings around each sphere, marked by crosses. This 
is the flow pattern expected if the flows of the two spheres 
do not interact strongly [29]. Under these conditions the 
time-averaged center of mass of the two spheres remains 
stationary: the dimer does not swim. 

As the amplitude increases, the flows grow in strength, 
but more importantly, the flows around each sphere start 
to interact strongly. The lower loop of the upper sphere 
is forced towards the surface of the sphere and reduces in 
size. Eventually, for sufficiently high amplitudes, there 
are only three vortex loops, as shown in Fig. 3 (c). A jet 
of fluid directed downwards from the smaller sphere can 
be observed from the plot of the normalised velocity field. 
Fig. 3 (c), and from experiment Fig. 3 (a). Under these 
conditions, the swimmer moves upwards, in the opposite 
direction to the jet. 

Simulations also allow us to vary parameters which 
are not easily accessible experimentally, such as a wider 
range of fluid viscosities, as shown in Fig. 4. When the 
dimer is moving there are four independent length scales: 
n//, A^, L and the viscous length 5. We obtain the best 


data collapse if n// is made non-dimensional by dividing 
by L rather than either of the other two length scales 
(see Supplemental Material [24]). Fig. 4 shows the simu¬ 
lation data plotted in this way indicating data collapse, 
the same way as the experimental data collapse shown in 
Fig. 1. The lower right inset to Fig. 4. shows typical 
trajectories after vibration has been applied. There are 
a few seconds of transient motion before the steady-state 
velocity is reached. 

Figs. 1, 2 and 4 all show that v/fL scales approx¬ 
imately linearly with the streaming Reynolds number 
Rcs for sufficiently large amplitudes A^. This behaviour 
is different from that observed for magnetic granular 
snakes m and rigid dimers on surfaces [31]. A sim¬ 
ple argument can be constructed to explain the scaling 
behaviour. Taking the unequal-sized swimmer as an ex¬ 
ample, the smaller sphere has a much larger amplitude 
of motion than the larger sphere, (see movie in the Sup¬ 
plemental Material [24]). The smaller sphere acts as a 
pump, imparting downward momentum to the fluid. The 
reaction force on the small sphere is equal and opposite 
to the rate of momentum transfer to the fluid. Its mag¬ 
nitude is proportional to the square of its speed (/A^)^ , 
the fluid density, p, and a geometric factor proportional 
to L^. In this simple model, the force is balanced by the 
Stokes’ drag on the larger sphere which scales as QttLtiv 
with V the velocity of the swimmer and rj is the dynamic 
viscosity of the fluid (pz^). By equating the two forces we 
obtain v/fL proportional to Rcs — as observed in 

the data for large amplitudes. 

Note, however, that this particular scaling behaviour 
is not expected to hold generally as there are four in¬ 
dependent length scales in this problem, and therefore 
three independent dimensionless ratios of lengths. The 
argument presented above is only expected to hold in the 
limit 5. 

The analysis given above assumes a strong asymme¬ 
try of the flows around both spheres, an assumption that 
breaks down at lower Reynolds numbers, as shown from 
the flow patterns in Fig. 3. In both experiment and 
simulation there appears to be a critical onset value of 
Res ~ 20 for swimming to occur, obtained by extrap¬ 
olation of the data to u = 0. It has been argued that 
asymmetric objects have a continuous transition to swim¬ 
ming m- This is not necessarily inconsistent with our 
observations. For Res below the apparent onset it is dif¬ 
ficult to determine whether v is strictly zero or is just 
small: experimentally it is hard to ensure that any small 
centre-of-mass motion is not due to residual buoyancy; 
in simulation, lattice effects may influence the motion 
when the amplitudes of the spheres become comparable 
to the lattice spacing. The existence of an apparent on¬ 
set to motion has also been observed in an asymmetric 
flapping wing m and the ‘acoustic scallop’ [20]. The 
good agreement between the experiment and simulation 
for our system allows us to conjecture that the apparent 
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onset of motion arises from the change in topology of the 
streaming flows. 

The examples presented here show a rich variety of 
behaviour but only represent a small part of the param¬ 
eter space. A systematic investigation into the influence 
of the overall size of the dimer, the ratio of the sphere 
diameters, the sphere density ratios and the gap width 
would be informative. It would be of interest to make a 
fully self-propelled swimmer based on the relative vibra¬ 
tion of two spheres, driven by an internal linear motor, 
since such swimmers would not be constrained to move 
along one axis. Collections of such swimmers might be 
expected to exhibit interesting cooperative behaviour in¬ 
duced by interacting streaming flows [30l [32l |33l [38l |39] . 
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